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Abstract 

Based on the distinction between the covariant and contravariant metric tensor 
components in the framework of the affine geometry approach and the s.c. "gravita- 
tional theories with covariant and contravariant connection and metrics" , it is shown 
that a wide variety of third, fourth, fifth, seventh, tenth- degree algebraic equations 
exists in gravity theory. This is important in view of finding new solutions of the 
Einstein's equations, if they are treated as algebraic ones. Since the obtained cubic 
algebraic equations are multivariable, the standard algebraic geometry approach 
for parametrization of two-dimensional cubic equations with the elliptic Weierstrass 
function cannot be applied. Nevertheless, for a previously considered cubic equa- 
tion for reparametrization invariance of the gravitational Lagrangian and on the 
base of a newly introduced notion of "embedded sequence of cubic algebraic equa- 
tions", it is demonstrated that in the multivariable case such a parametrization is 
also possible, but with complicated irrational and non-elliptic functions. After find- 
ing the solutions of a system of first - order nonlinear differential equations, these 
parametrization functions can be considered also as uniformization ones (depending 
only on the complex uniformization variable z) for the initial multivariable cubic 
equation. 



1 INTRODUCTION 



Inhomogeneous cosmological models have been intensively studied in the past in reference 
to colliding gravitational waves [1] or singularity structure and generalizations of the Bondi 
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- Tolman and Eardley-Liang-Sachs metrics [2, 3]. In these models the inhomogeneous 
metric is assumed to be of the form [2] 

ds'' = - e2"(*''-'f'^)dr2 - e2'^(*''-'^''^)(ci?/^ + dz'') (1.1) 

(or with r ^ z and z ^ x), which is called the Szafron-Szckeres metric [4-7]. In [7], after 
an integration of one of the components - Gi of the Einstein's equations, a solution in 
terms of an elliptic function is obtained. 

In different notations, but again in the framework of the Szafron-Szekerez approach 
the same integrated in [7] nonlinear differential equation 

^ + 2M(^)$-^ + U^' (1.2) 

was obtained in the paper [8] of Kraniotis and Whitehouse. They make the useful obser- 
vation that (1.2) in fact defines a (cubic) algebraic equation for an elliptic curve, which 
according to the standard algebraic geometry prescribtions (see [9] for an elementary, 
but comprehensive and contemporary introduction) can be parametrized with the elliptic 
Weierstrass function 

/ N 1 V- r 1 1 " 



[z — ujS^ uS^ 



(1.3) 



and the summation is over the poles in the complex plane. Two important problems 
immediately arise, which so far have remained without an answer: 

1. The parametrization procedure with the elliptic Weierstrass function in algebraic 
geometry is adjusted for cubic algebraic equations with number coefficients! Unfortu- 
nately, equation (1.2) is not of this type, since it has coefficient functions in front of the 
variable $, which depend on the complex variable z. In view of this, it makes no sense 
to define "Weierstrass invariants" as 

y2 12 ' 216 12 ^ > ' ^ > 

since the above functions have to be set up equal to the complex numbers gi and g-^ (the 
s. c. Eisenstein series) 

^2 = 60V4 = y"7 , (1-5) 

^3 = i4oy^ = y"7 — ^- — (1.6) 

^•^ Z^ Z^ (n + mT)6 ^ ' 

u)<ZX n,m ^ ' 

and therefore additional equations have to be satisfied in order to ensure the parametriza- 
tion with the Weierstrass function. 

2. Is the Szekerez - Szafron metric the only case, when the parametrization with the 
Weierstrass function is possible? Closely related to this problem is the following one - is 
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only one of the components of the Einstein's equation parametrizable with p{z) and its 
derivative? 

This paper has the aim to present an adequate mathematical algorithm for finding 
solutions of the Einstein's equations in terms of elliptic functions. This approach is based 
on the clear distinction between covariant and contravariant metric tensor components 
within the s.c affine geometry approach, which will be clarified further in Section 2. Af- 
terwords, a cubic algebraic equation in terms of the contravariant metric components 
will be obtained, which according to the general prescription and the algorithm in the 
previous paper [10] can be parametrized with the Weierstrass function and its derivative. 
Respectively, if the contravariant components are assumed to be known, then a cubic 
(or a quartic) algebraic equation with respect to the covariant components can be in- 
vestigated and parametrized again with the Weierstrass function. Thus it will turn out 
that the parametrization with the Weierstrass function will be possible not only in the 
Szafron-Szekeres case, but also in the general case due to the "cubic" algebraic structure 
of the gravitational Lagrangian. This is an important point since valuable cosmological 

characteristics for observational cosmology such as the Bubble's constant H{t) = and 

the deceleration parameter q = j^-^^y expressed in terms of the Jacobi's theta 

R it) 

function and of the Weierstrass elliptic function respectively [8]. Unfortunately, in the 
paper [8] the Eisenstein series (1.5-1.6) have not been taken into account, due to which 
the obtained expression for the metric will be another one and will be modified. 

Instead of searching out eliptic solutions of the Einstein's equations for each separate 
case of a given metric, as in nearly all of the mentioned papers, in the this paper an- 
other method will be proposed. First, a cubic algebraic equataion will be parametrized 
with respect to one of the contravariant components, following the approach in a pre- 
vious paper [10]. Then, this parametrization will be extended to more than one vari- 
able in the multivariable cubic algebraic equation (section 6). This will be a substan- 
tial and new development, different from the standard algebraic geometry approach, in 
which only two-dimensional cubic equations are parametrized with the (elliptic) Weier- 
strass function and its derivative. Finally, the dependence of the generalized coordinates 

— X^{xi,X2,X3, ,Xn) on the complex variable z will be established from a derived 

system of first-order nonlinear differential equations (section 7). The generalized coor- 
dinates can be regarded as n— dimensional hypersurfaces, defining a transition from an 
initiallv defined set of coordinates chosen manifold to another set 

of the generalized coordinates ,X^. Since the covariant metric components 

gij also depend on these coordinates, this means that their dependence on the complex 
variable z will also be known. In other words, at the end of the applied approach, each 
initially given function gij{t,x^) of the time and space coordinates will be expressed also 
as gij{z). The algebraic approach will be applied to the s .c. cubic algebraic equation 
for reparametrization invariance of the gravitational Lagrangian, but further it will be 
shown that not only the approach will be applicable in the general case of an arbitrary 
contravariant tensor, but also concrete solutions for the metric gij{z) will be given in the 
case of specially chosen simple metrics. 
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The present paper continues and develops further the approach from a previous paper 
[10], where a definite choice of the contravariant metric tensor was made in the form of the 
factorized product g"^^ = dX^dXK The differentials dX^ are assumed to lie in the tangent 
space Tx of the generalized coordinates. In Section 2 of the present paper some basic 
facts about the affine geometry approach and the s.c. gravitational theory with covariant 
and contravariant metrics and connections ( GTCCMC) will be reminded, which has been 
described in the review article [13]. In its essence, the distinction between covariant and 
contravariant components is related to the affine geometry approach [15, 16], according 
to which the four-velocity tangent vector at each point of the observer's worldline is not 
normalized and not equal to one, i.e. IJ,"" — P ^ 1. Similarly, for a second-rank tensor 
one would have g^^g'^"' — ^'jl- ^^e next section 3 it will be demonstrated briefly 
how the cubic algebraic equation with respect to the differentials dX^ was derived in 
[10], but in fact the aim will be to show that depending on the choice of variables in 
the gravitational Lagrangian or in the Einstein's equations, a wide variety of algebraic 
equations (of third, fourth, fifth, seventh degree) in gravity theory may be treated, if a 
distinction between the covariant metric tensor components and the contravariant ones is 
made. This idea, originally set up by Schouten and Schmutzer, was further developed 
in the papers [13, 14]. In usual gravity theory, the contravariant components are at the 
same time inverse to the covariant ones , and thus the correspondence between " covectors" 
(in our terminology - these are the "vectors") and the "vectors" (i.e. the contravariant 
vectors") is being set up, since both these kinds of tensors satisfy the matrix equation 
9ii9^^ = ^i^- However, within the framework of affine geometry, such a correspondence is 
not necessarily to be established (see again [15-18]) and both tensors have to be treated 
as different mathematical objects, defined on one and the same manifold. 

The physical idea, which will be exploited in this paper will be: can such a gravitational 
theory with a more general contravariant tensor have the same gravitational Lagrangian 
as in the known gravitational theory with contravariant metric tensor components, which 
are at the same time the inverse ones to the covariant one? On the base of such an 
"equivalence" the s. c. cubic algebraic equation for reparametrization invariance of the 
gravitational Lagrangian was obtained in [10]. The derivation was based also on the 
construction of another connection = \dX^dX^{gik,i+gii,k~9ki,i)- It can be proved that 
the connection F^^ has two very useful properties: 1. It may have an affine transformation 
law under a broad variety of coordinate transformations, which can be found after solving 
a system of nonlinear differential equations. 2. F^^ is an equiaffine connection, which 
is a typical notion, introduced in classical affine geometry [15, 16] and meaning that 
there exists a volume element, which is preserved under a parallel displacement of a basic 
n— dimensional vector e = ej^jj. -jn- Equivalently defined, F^^ is an equiaffine connection 
[15, 16] if it can be represented in the form F^^ = dklge, where e is a scalar quantity. This 
notion turns out to be very convenient and important, since for such types of connections 
we can use the known formulae for the Ricci tensor, but with the connection F^^ instead 
of the usual Christoffell one F|^. Moreover, the Ricci tensor R^j will again be a symmetric 
one, i.e. R., = i?,, = - a.F^ ■ + fl,f\^ - F-f 
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2 AFFINE GEOMETRY APPROACH AND GRAV- 
ITATIONAL THEORIES WITH COVARIANT 
AND CONTRAVARIANT METRICS AND CON- 
NECTIONS 

This section has the purpose to review some of the basic aspects of GTCCMC, which 
would further allow the application of algebraic geometry and theory of algebraic equations 
in gravity theory. 

It is known in gravity theory that the metric tensor g^j determines the space - time 
geometry, which means that the Christoffell connection 

r'fe = i^g^^{gks,i + gis,k - 9ik,s) (2.1) 

and the Ricci tensor 

D ik il I -pi -pm -pm-pl (n r\\ 

can be calculated. 

It is useful to remember also [23] the s. c. Christoffell connection of the first kind: 

^i;kl = ginJ^H = -;^{gik,l + gil,k — gkl,i) , (2.3) 

obtained from the expression for the zero covariant derivative = Vigik — gik,i — 5'm(«r^; 
. By contraction of (2.3) with another contravariant tensor field g^^, one might as well 
define another connection: 

^ki = g^^^i-M = g^^gimTTi = -jg^^{gik,i + gii,k — gki,i) , (2.4) 

not consistent with the initial metric gij. Clearly the connection (2.4) is defined under the 
assumption that the contravariant metric tensor components g^^ are not to be considered 
to be the inverse ones to the covariant components gij and therefore g^^gim = /m(^)- 

In fact, the definition g^^gim = turns out to be inherent to gravitational physics. 
For example, in the projective formalism one decomposes the standardly defined metric 
tensor (with gijg^'' — S'^) as 

gij = Pij + hij , (2.5) 

together with the additional assumption that the two subspaces, on which the projective 
tensor pij and the tensor hij are defined, are orthogonal. This means that 

Pijh^'' = . (2.6) 

As a consequence 

Pijpi'^ ^ 51 - Kjh^' ^ 51 , (2.7) 
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meaning that the contravariant projective metric components jp^ (in the orthogonal sub- 
space to the tensor hij) are no longer inverse to the covariant ones pij. 

An example of gravitational theories with more than one connection are the so called 
theories with affine connections and metrics [13], in which there is one connection V^^^ for 
the case of a parallel transport of covariant basic vectors Vg^Ca = T^^p and a separate 
connection Pj^ for the contravariant basic vector e''', the defining equation for which is 
Ve^e" = P"^ . In these theories, the contravariant vector and tensor fields are assumed 
to be not the inverse ones to the covariant vector and tensor fields. This implies that 

ey ^ /f (x) ^ 5i (2.8) 

and consequently, a distinction is made between covariant and contravariant metric ten- 
sors (and vectors too). Clearly, in the above given case (2.7) of projective gravity, this 
theory should be considered as a GTCCMC. In the same spirit, since the well - known 
Arnowitt - Deser - Misner (ADM) (3-1-1) decomposition of spacetime [43, 44] is built upon 
the projective transformation (2.5), it might be thought that it should also be considered 
as such a theory. But in fact, the ADM (3+1) formalism definitely is not an example for 
this due to the special identification of the vector field's components [43, 44] w^^h, certain 
components of the projective tensor, as a consequence of which all the (spacelike defined) 
contravariant projective tensor components p"'^ (a,j3,^ — 0,1,2,3 ; i,j = 1,2,3^ turn 
out to be the inverse ones to the covariant projective components Pa-y- In the case of the 
ADM (3-1-1) decomposition, such an identification is indeed possible and justified, since 
the gravitational field posseses coordinate invariance, allowing to disentangle the dynam- 
ical degrees of freedom from the gauge ones. But in the case when the tensors hij are 
related with some moving matter (with a prescribed motion) and an observer, "attached" 
to this matter "measures" all the gravitational phenomena in his reference system by 
means of the projective metric pij, this will be no longer possible. Then the relation 
(2.7) will hold, and the resulting theory will be a GTCCMC. Naturally, if the tensor hij 
in (2.5) and (2.7) is taken in the form hij = -UiU^ and if the vector field u (tangent 
at each point of the trajectory of the moving matter) is assumed to be non- normalized 
(i.e. e{x) = UiU^ ^ 1), then one would have to work not within the standard relativistic 
hydrodynamics theory (where pij = gij — UiUj and Pijp'^ = ~ UiU^), but within the 
formalism of GTCCMC (where Pijp''' = fi = 5f — ^Uiu'^ ^ S'l). One may wonder why 
this should be so, since the last two formulaes for Pijp''^ for the both cases look very much 
alike, with the exception of the "normalization" function ^ in the second formulae. But it 
shall become clear that in the first case the right-hand side has a tensor transform,ation 
property, while in the second case due to the function ^ there would be no such property. 
And this shall turn out to be crucial. 

In order to understand this also from another point of view, let us perform a covariant 
differentiation of both sides of the relation (2.8). Then one can obtain that the two 
connections are related in the following way [13] 

flk = ^'jkfi + pui ; ifik = dj;) . (2.13) 
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Note also the following important moment - (x) are considered to be the components 
of a function. Otherwise, if they are considered to be a (mixed) tensor quantity, the 
covariant differentiation of the mixed tensor f^{x) in the right-hand side of e^e'^ = f^{x) 
would give exactly the same quantity as the one in the left-hand side for every choice of 
the two connections F^^ and P/^, including also for the standard case (Einsteinian gravity) 

= ~r'fc. This would mean that from a mathematical point of view there would be no 
justification for the introduction of the second covariant connection P/^ . However, since 
fa{x) are related with the description of some moving matter in the Universe, a tensor 
transformation law should not be prescribed to them. So they should remain components 
of a function and consequently, the introduction of the second connection P^*^ is inevitable. 

In order to understand further why and in what cases the distinction between covariant 
and contravariant metric components will lead to an inevitable introduction of two different 
connections F^j and Plj, let us prove the following statement: 

Proposition 1 // ci, 62, e„ is a basis of covariant vector fields and are the com- 
ponents of a given function or a constant, then a basis of contravariant basic fields 
~ai —02^ —an gQj^ found so that for each i and aj one has e^e"^ = . 

This proposition is in fact is a generalization of the well-known theorem from differ- 
ential geometry that if a basis of (covariant) vector fields is given, then a dual basis of 
(contravariant) vector fields can be found, so that the contravariant vector fields are the 
inverse ones to the covariant ones, i.e. Cje"-' = 5". 

The proof is very simple, but essentially based on the relation (2.13). If the covariant 
basic vector fields are given, then the contravariant connection components F^^ will be 
known too. Since /j^, are derivatives of a function, one may take the expression (2.13) 

/j J. = r';,, fi+Pii^fj , which for the moment shall be treated as a system of linear 
algebraic equations with respect to the (unknown) connection components Pi\. A solution 
of this system can be found for the connection components P;*^. Then the condition for the 
parallel transport of the contravariant basic vector fields Ve^e" = P"^ e^ can be written 
as dpe"' — and considered as a system of n ordinary differential equations with 

respect to the components e". From the solution of this system, e'" can be found. So the 
new basis e"^, e"" will be uniquely determined, up to the integration constants, 
contained in the solution of the system of differential equations. 

After proving this proposition, the difference between standard relativistic hydrody- 
namics and "modified" relativistic hydrodynamics with a variable length can be easily 
understood. In the first case, the right-hand side in pijP^^ = ~ ""i^'^ = ft 7^ trans- 
forms as a tensor, which is ensured also by normalization property -Uj-u* = 1. Therefore 
(2.13) and the proposition will not hold, so the contravariant basic vector fields are de- 
termined in the standard way e^e^ — and more importantly, they cannot be determined 
in another way, in spite of the fact that fi ^ 

In the second case, the situation is just the opposite - the right-hand side of Pijp'^ = 
~ = fi 7^ transforms not as a tensor because of the "normalization" factor 
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^, the proposition holds and thus the basic vector fields are determined as eie^ = //. 
Therefore, the treatment of relativistic hydrodynamics with "variable length" should be 
withm the GTCCMC. 

In the present case, the introduced new connection (2.4) should not he identified with 
the connection Pj^, since the connection V^i = g'^^T^-ki is introduced by means of mod- 
ifying the contravariant tensor and not on the base of any separately defined parallel 
transport for the contravariant basic vectors. Moreover, the connection T'li turns out to 
be a linear combination of the Christoffell connection components F^^, and the relation 
between them is not of the type (2.13). In such a way, there will not be a contradiction 
with the case when the two connections F^^ and F^^ are not defined as separate ones, 
since later on, in deriving the cubic algebraic equation in the general case and for the 
case g^'' = dX^dX^ also, it would be supposed that 5^** is a tensor. This would mean 
(from g^'^gim = fmi^)) ^^^^ /m(^) ^^^^ ^^^^ t)e a (mixed) tensor quantity, and therefore 
the covariant differentiation of e^e''^ = /^(a;) will not produce any new relation. 

3 BASIC ALGEBRAIC EQUATIONS IN GRAV- 
ITY THEORY. TENSOR LENGTH SCALE 

If one writes down the Ricci tensor in terms of the newly defined contravariant tensor 
g^^ = dX^dX^ , then the following fourth - degree algebraic equation can be obtained 



Rik — dX 



d{dX') _ 1 1 d{dX') 



+ 



+ -dX''dX"^dX^dX^ [gm[k,t9l]r,i + gi[l,tgmr,k] + '^9t[k,i9mr,l]\ , (3.1) 

where p is the scalar quantity 

p = div{dX) = ^^, (3.2) 

which "measures" the divergency of the vector field dX. The algebraic variety of the 
equation consists of the differentials dX^ and their derivatives ^^^§p- ■ 

In the same spirit, one can investigate the problem whether the gravitational La- 
grangian in terms of the new contravariant tensor can be equal to the standard repre- 
sentation of the gravitational Lagrangian. This standard (first) representation of the 
gravitational Lagrangian is based on the standard Christoffell connection F^^- (given by 
formulae (2.1)), the Ricci tensor Rtk (formulae (2.2)) and the other contravariant tensor 
g'^ = dX'dX^ [10] 

Li^-sf^g"'Rik^-V^dX'dX^Rik . (3.3) 

In the second representation, the Christoffell connection F^^ and the Ricci tensor R^^ 
are "tilda" quantities, meaning that the "tilda" Christoffell connection is determined by 
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formulae (2.4) with the new contravariant tensor g'^^ — dX'^dX^ and the "tilda" Ricci 
tensor Rik - by formulae (2.2), but with the "tilda" connection F^^- instead of the usual 
Christoffell connection Tfj. Thus the expression for the second representation of the 
gravitational Lagrangian acquires the form 

L2 = -V^g'^Rii = -V^dX'dX'{prigkrdX'' - ri^gi^d^X' - ^^.g^yd'^X'^} . (3.4) 

The condition for the equivalence of the two representations Li = L2 gives a cubic alge- 
braic equation with respect to the algebraic variety of the first differential dX'^ and the 
second ones di^X^ [10] 

dX'dX' {pTligkrdX'^ - Tl.gird^X'^ - Tl^.g^yd^X^) - dX'dX'Ra = . (3.5) 

In [22] also another cubic algebraic equation has been obtained, but after the application 
of a variational approach. 

Following the approach in [10], the Einstein's equations in vacuum for the general 
case were derived under the assumption that the contravariant metric tensor components 
are the "tilda" ones: 

~ 1 ~ ~ 1 m n~ 

= Rij — -gijR — Rij — 29ijdX dX Rmn — 

= -lpgijr:^n9krdX'dX"^dX- + \gij{Tl^gnr + K^m9k)r)d''XHX^dX^+ 

+pV\^gkrdX^ - {Tl^g^r + V]^^^gk)r)d^X^ . (3.6) 

This equation represents again a system of cubic equations. In addition, if the differ- 
entials dX* and d'^X'^ are known, but not the covariant tensor gij^ the same equation can 
be considered also as a cubic algebraic equation with respect to the algebraic variety of 
the metric tensor components gij and their first derivatives gij^k- 

It might be thought that the definite choice of the contravariant tensor in the form of 
the factorized product 'g'^^ — dX^dX^ is a serious restriction, in view of the fact that the 
second derivatives of the covariant tensor components gij^ki are not present in the equation. 
This is indeed so, because the algebraic structure of the equation is simpler to deal with 
in comparison with the general case, and so it is easier to implement the algorithm for 
parametrization, developed in [10]. But there is one argument in favour of this choice 
(although the case for an arbitrary contravariant tensor is no doubt more important) - 
since the metric can be expressed as ds^ = l{x) = QijdX^dX^ (consequently dX^dX^ = 
l{x)g^^), the obtained cubic algebraic equations (3.5) and (3.6) can be considered in regard 
also to the length function l{x). Since for Einsteinian gravity gijg^'^ = S'^ (i.e. g^'^ = 'g^^ — 
dX^dX^), then for this case the length function is "postulated" to be Z = 1. But the 
length function can also be obtained as a solution of the cubic equation, and thus in 
more general theories of gravity solutions with I 7^ 1 may exit. In fact, for a general 
contravariant tensor 'g'^^ ^ dX^dX^, one would have g'^^ — llg'^-' , where will be proposed 
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to be called a "tensor length scale", and the previously defined length function l(x) is a 
partial case of the tensor length scale for /* = The physical m,eaning of the notion of 
tensor length scale is simple - in the different directions (i.e. for different i and j) the length 
scale is different. In particular, some motivation for this comes from Witten's paper [45], 
where in discussing some aspects of weakly coupled heterotic string theory (when there 
is just one string couplings ) and the obtained too large bound on the Newton's constant 
it was remarked that "the problem might be ameliorated by considering an anisotropic 
Calabi - Yau with a scale \f(^ in d directions and j^^^ in (6 — d) directions". Thus it 
may be proposed to realize this if one takes 

= gijdXUx'' = for i, j. A; = 1, d , (3.7) 

ll = gacdX^dX^ = 128^ ioY a,h,c = d + l,....,Q . (3.8) 

Note also the justification for the name "tensor length scale" - if l\ is a tensor quantity, 
so will be the "modified" contravariant tensor 'g'^^ = l\g^^ , and consequently in accord 
with section 2 there will be no need for the introduction of a new covariant connection 

Pp.. And this is indeed the case, because the relation between the two connections T^, and 

''J ^ 

is given by formulae (2.4) F^^ := 'g'^girrX^i- In other words, these two connections are 
not considered to be "separately introduced" and so they do not depend on one another 
by means of the equality (2.13). 

The purpose of the present paper further will be: how can one extend the proposed in 
[10] approach for the "modified" contravariant metric components (as g'^^ — dX^dX^) to 
the case of a generally defined contravariant tensor g^^ ^ dX^dX^ ? 

4 INTERSECTING ALGEBRAIC VARIETIES AND 
STANDARD (EINSTEINIAN) GRAVITY THE- 
ORY 

A more general theory with the definition of the contravariant tensor as g^^ = dX^dX^ 
should contain in itself the standard gravitational theory with gijg^^ = 6^. From a 
mathematical point of view, this should be performed by considering the intersection [19, 
20, 21] of the cubic algebraic equations (3.6) with the system of quadratic algebraic 
equations for the algebraic variety of the n variables 

gijdXUX^ = 5f . (4.1) 

In its general form gifg^^ — 5\ with an arbitrary contravariant tensor 'g^^, this system can 
also be considered together with the Einstein's "algebraic" system of equations, which in 
the next section shall be derived for a generally defined contravariant tensor. From an alge- 
bro - geometric point of view, this is the problem about the intersection of the Einstein's 

algebraic equations with the system of (linear) hypersurfaces for the ( ^ | + n 
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contravariant variables, if the covariant tensor components are given. Since the derived 
Einstein's algebraic equations are again cubic ones with respect to the contravariant met- 
ric components, this is an analogue to the well - known problem in algebraic geometry 
about the intersection of a (two-dimensional) cubic curve with a straight line. However, 
in the present case the straight line and the cubic curve are multi - dimensional ones, 
which is a substantial difference from the standard case in algebraic geometry. 

The standardly known solutions of the Einstein's equations can be obtained as an 
intersection variety of the Einstein's algebraic equations with the system Qijcj^^ = 6^ . 
However, the strict mathematical proof that such an intersection will give the known 
solutions is still lacking. 

5 ALGEBRAIC EQUATIONS FOR A GEN- 
ERAL CONTRAVARIANT METRIC TENSOR 

Let us write down the algebraic equations for all admissable parametrizations of the 
gravitational Lagrangian for the generally defined contravariant tensor g^^ , following the 
same prescription as in section 3, where the equality of the two representations of the 
gravitational Lagrangian has been supposed: 

+rrr'9pr9,s (n^r^ - rlru -R^o . (5.1) 

This equation is again a cubic algebraic equation with respect to the algebraic variety of 
the variables g^^ and g^l, and the number of variables in the present case is much greater 
than in the previous case for the contravariant tensor g^^ = dX'^dX^ . At the same time, 
this equation is a fourth - degree algebraic equation with respect to the covariant metric 
tensor gij and its first and second partial derivatives. With respect to the algebraic variety 
of all the variables g'\ g'l, gij, gij^ki 9ij,kh the above algebraic equation is of seventh order 
and with coefficient functions, due to the presence of the terms with the affine connection 
r?^ and its derivatives, which contain the contravariant tensor g''^ and g^l- 

Similarly, the Einstein's equations can be written as a system of third - degree alge- 
braic equations with respect to the (generally chosen) contravariant variables and their 
derivatives ^ 

= Rij — -gijR = 

~9 ^r;i\j)-il]'^9 \]yr;ij]'^ 9 9 iX r;i3^ s;lm ~ ^ s;il^ r;krn) ~ 
-\9^r9'^^'tlTl,,grs - \9^r9'^'''9'^' {^'^k9rs),i - 

-\9iir'r9'^'9^9,s{Kk'^L-Kinj ■ (5-2) 

Interestingly, the same system of equations can be considered as a system of fifth - degree 
equations with respect to the covariant variables (which is the difference from the previous 



11 



case). The mathematical treatment of fifth - degree equations is known since the time 
of Fehx Klein's famous monograph [24], published in 1884. A way for resolution of such 
equations on the base of earlier developed approaches by means of reducing the fifth - 
degree equations to the so called modular equation has been presented in the more recent 
monograph of Prasolov and Solov'yev [9]. Some other methods for solution of third-, 
fifth- and higher- order algebraic equations have been given in [25, 26]. A complete 
description of elliptic, theta and modular functions has been given in the old monographs 
[27, 28]. Also, solutions of n— th degree algebraic equations in theta - constants [29] and 
in special functions [30] are interesting in view of the not yet proven hypothesis in the 
paper by Kraniotis and Whitehouse [8] that "all nonlinear solutions of general relativity 
are expresed in terms of theta - functions, associated with Riemann - surfaces". Some 
other monographs, related to elliptic functions and elliptic curves are [31-42]. 
Two other important problems can be pointed out: 

1. One can find solutions of the system of Einstein's equations not as solutions of a 
system of nonlinear differential equations, but as elements of an algebraic variety, satisfy- 
ing the Einstein's algebraic equations. The important new moment is that this gives an 
opportunity to find solutions of the Einstein's equations both for the components of the 
covariant metric tensor gij and for the contravariant ones 'g^^. This means that solutions 
may exist for which gijg^^ 7^ • In other words, a classification of the solutions of the 
Einstein's equations can be performed in an entirely new and nontrivial manner - under a 
given contravariant tensor, the covariant tensor and its derivatives have to be found from 
the algebraic equation, or under a given covariant tensor, the contravariant tensor and its 
derivatives can be found. 

2. The condition for the zero - covariant derivative of the covariant metric tensor 
^kQij = and of the contravariant metric tensor Vug''-' = can be written in the form of 
the following cubic algebraic equations with respect to the variables gij, gij^k and 'g^^ : 

^kQij = gij,k - ^k{i9j)i = 9ij,k - g^''^s;k{i93)i = (5.3) 

and 

Q^Vkf^ ^fi + J^Y^'rr.,sk . (5.4) 

The first equation (5.3) is linear with respect to g^'* and quadratic with respect to g^, 
gij^k, while the second equation (5.3) is linear with respect to g^j, gij^k a-nd quadratic with 
respect to g^^. 

6 EMBEDDED SEQUENCE OF ALGEBRAIC 
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EQUATIONS AND FINDING THE SOLU- 
TIONS OF THE CUBIC ALGEBRAIC EQUA- 
TION 

The purpose of the present subsection will be to describe the method for finding the 
solution (i. e . the algebraic variety of the differentials rfX*) of the cubic algebraic equation 
(3.5) (in the limit d'^X'' = 0). The method has been proposed first in [10] but here it will 
be developed further and applied with respect to a sequence of algebraic equations with 
algebraic varieties, which are embedded into the initial one. This means that if at first 
the algorithm is applied with respect to the three-dimensional cubic algebraic equation 
(3.5) and a solution for rfX^ (depending on the Weierstrass function and its derivative) is 
found, then the same algorithm will be applied with respect to the two-dimensional cubic 
algebraic equation with variables dX^ and dX^, and finally to the one-dimensional cubic 
algebraic equation of the variable dX^ only. 

The basic and very simple idea about paramctrization of a cubic algebraic equation 
with the Weierstrass function [9, 11,12] can be presented as follows: Let us define the 
lattice A = {mui + nu2 \ m,n E Z; Ui,uj2 E C, Im^ > 0} and the mapping / : C/A ^ 
CP^, which maps the factorized (along the points of the lattice A) part of the points on 
the complex plane into the two dimensional complex projective space CP^. This means 
that each point z on the complex plane is mapped onto the point {x,y) — {p{z),p (z)), 
where x and y belong to the affine curve 

= 4:x^ - g2x - 93 ■ (6.1) 

In other words, the functions x = p{z) and y = p (z), where p{z) denotes the Weierstrass 
elliptic function (1.3), are uniformization functions for the cubic curve. It can be proved 
[9] that the only cubic algebraic curve with number coefficients, which is parametrized by 
the uniformization functions x — p{z) and y — p (z) is the affine curve (6.1). 

In the case of the cubic equation of reparametrization invariance (3.5), the aim will be 
again to bring the equation to the form (6.1) and afterwards to make equal each of the 
coefficient functions to the (numerical) coefficients in (6.1). 

In order to provide a more clear description of the developed method, let us divide it 
into several steps. 

Step 1. The initial cubic algebraic equation (3.5) is presented as a cubic equation 
with respect to the variable dX^ only 

A^idX')' + B-i{dXY + C-iidX') + G^^HdX^, dX\ gi,, TJ,-, Rh,) = , (6.2) 

where naturally the coefficient functions A3, B-^ , C3 and G^"^^ depend on the variables 
dX^ and dX'^ of the algebraic sub variety and on the metric tensor gij, the Christoffel 
connection T^j and the Ricci tensor Rij: 

As = 2pri^g3r ; S3 = 6pr;3^3,dX- - i?33 , (6.3) 
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Ca = -2R^^dX'' + 2p(r;^^3r + 2Tlpg^,)dX''dXP . (6.4) 

The Greek indices a, (3 take values a, (3 — 1,2 while the indice r takes values r = 1, 2, 3. 
Step 2. A linear- fractional transformation 



3 _ a3,{z)dX + 63(2;) 



cz{z)dX +d^{z) 



(6.5) 



is performed with the purpose of setting up to zero the coefficient functions in front of 

3 

the highest (third) degree of dX . This will be achieved if G^'^\dX'^, dX'^, g^j, V^^ Rik) = 
where 

Q = Asal + C3CI + B^a^cs + 2c^d^C^ . (6.6) 

This gives a cubic algebraic equation with respect to the two-dimensional algebraic variety 
of the variables dX^ and dX'^: 

pVl^^^gpydX^dX'^dX^ + K^^ldX'^dX^ + + 2p Vl^g^^ = (6.7) 

and K^^l and K^a^ are the corresponding quantities [10] 



= + 2p-(i + 2-)(2r;^^3. + n^g^r 
C3 C3 



(6.8) 



and 



C3 



3p-r;353r-(l + 2^)i?«3 
C3 C3 



(6.9) 



Note that since the linear fractional transformation (with another coefficient functions) 

will again be applied with respect to another cubic equations, everywhere in (6.5 - 6.8) the 

coefficient functions 03(2;), 63 (-2), c^{z) and d^{z) bear the indice "3", to distinguish them 

from the indices in the other linear-fractional tranformations, which are to be applied. In 

— 3 

terms of the new variable — dX the original cubic equation (3.5) acquires the form 
[10] 



= Pi(n) mr + P2{n) + P^{n) m + Pi{n) , 



(6.10) 



where Pi(n) ,P2{n), P^iji) and P4(n) are comphcated functions of the ratios ^, ^ and 
A3, S3, C3 (but not of the ratio ^, which is very important). The variable m denotes the 
ratio ^ and the variable n is related to the variable n through the expresssion 



n — 



(^3 



5, 



(3) 



(6.11) 



where 



; _ ^3 C3 , C3 

dsds ds 



(6.12) 
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The subscript "3" in L^^ and L^^^ means that the corresponding ratios in the right- 
hand side also have the same subscript. Setting up the coefficient functions Pi{n) , P2(?^), 
P-i{n) equal to the number coefficients 4, 0, —g2, —Qs respectively, one can now parametrize 
the resulting equation 

= 4m^ — — (6.14) 
according to the standard prescription 

n^p{z) = ^ ; —=m^p{z) . (6.15) 

az C3 

Taking this into account, representing the linear-fractional transformation (6.5) as (divid- 
ing by C3) 

dX^ = ^ (6.16) 

dX +^ 

and combining expressions (6.11) for n and (6.16), one can obtain the final formulae for 
dX^ as a solution of the cubic algebraic equation 



- I P(^) r(3)B3 r(3) " ^^"^^^ 



In order to be more precise, it should be mentioned that the identification of the functions 
Pi(n) ,P2{n), Psin) with the number coefficients gives some additional equations [10], 
which in principle have to be taken into account in the solution for dX^. This has been 
investigated to a certain extent in [10], and will be continued to be investigated. Here 
in this paper the main objective will be to show the dependence of the solutions on the 
Weierstrass function and its derivative. Since only the ratios ^ and ^ enter these additional 
relations, and not ^ (which is related to the Weierstrass function), they do not affect the 
solution with respect to p{z) and p'{z). 

Since ^3 and C3 depend on dX^ and dX'^, the solution (6.17) for dX^ shall be called 
the embedding solution for dX^ and dX^. 

Step 3. Let us now consider the two-dimensional cubic equation (6.6). Following the 
same approach and finding the "reduced" cubic algebraic equation for dX^ only, it shall 
be proved that the solution for dX"^ is the embedding solution for dX^. 

For the purpose, let us again write down eq. (6.6) in the form (6.2), singling out the 
variable dX^: 

A2{dXY + B2{dXY + C2{dX^) + G^^\dX\gij,V\^,Rik) = , (6.18) 
where the coefficient functions ^42, i^2, C2 and G*^^) are the following: 

A2 = 2pVl^g2r ; B2 = K^^ +2p[2V\^g2r + Tl29iMX^ , (6-19) 
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C2 = 2p[r\,g2r + 2^{^g^r){dX'T + {K^ + K^^)dX^ + K?^ , (6.20) 

= 2pT[-,giridXy + Ki'^idX^y + xfUx^ + 2'pp\z)Vl^g^r • (6.21) 

Note that the starting equation (6.7) has the same structure of the first terms, if one 
makes the formal substitution —Raji K'^al second terms, but eq. (6.7) has 

two more additional terms K'^^dX^ + 2pp^[z)V'^2,9^r- Therefore, one might guess how 
the coefficient functions will look like just by taking into account the above substitution 
and the contributions from the additional terms. Revealing the general structure of the 
coefficient functions might be particularly useful in higher dimensions, when one would 
have a "chain" of cubic algebraic equations. Concretely for the three-dimensional case, 
investigated here, C2 in (6.20) can be obtained from C3 in (6.4), observing that there will 
be an additional contribution from the term Ka dX°' for a = 2. Also, in writing down 
the coefficient functions in (6.2) it has been accounted that as a result of the previous 
parametrization ^ = p{z) . 

Since eq. (6.18) is of the same kind as eq. (6.2), for which we already wrote down the 
solution, the expression for dX"^ will be of the same kind as in formulae (6.17), but with 
the corresponding functions A2,B2, C2 instead of ^43, S3, C3. Taking into account (6.19 - 
6.20), the solution for dX'^ can be written as follows: 

^^2 ^ ^P(^)P (^)V^+ /^l(^X^)^ + h2idX^) + 

^/{z)VC~2 + h{dX^)^ + l2{dX^) + h ' ^' ^ 

where hi, ^2, ^3, ^i, h, h are expressions, depending on ^, ^, F^^ (r = 1, 2, 3 ; a, /? = 1, 2), 

gap, , K2^^ and on the Weierstrass function. 

The representation of the solution for dX"^ in the form (6.22) shows that it is an 
"embedding" solution of dX^ in the sense that it depends on this function. Correspondingly, 
the solution (6.17) for dX^ is an "embedding" one for the variables dX^ and dX'^. 

Step 4. It remains now to investigate the one-dimensional cubic algebraic equation 

AiidX^Y + BiidXy + CiidX') + G^^\gi^, V%, R^k) = , (6.23) 

obtained from the two-dimensional cubic algebraic equation (6.18) after applying the 
linear-fractional transformation 

dX^ = ^-i^)dx\b2{z) ^g_24) 
C2{z)dX +d2{z) 

and setting up to zero the coefficient function before the highest (third) degree of {dX'^)^. 
Taking into account that as a result of the previous parametrization — = p{z) , the 
coefficient functions ^i,Si,Ciand Di are given in a form, not depending on dX'^ and 
dX^: 

Ai = 2pT[,gir , (6.25) 
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B, = Fsp{z) + Ki'^ = 2p(l + 2^) [2r[,g,r + ri,g,r]p{z) + k['^ , (6.26) 

+(1 + 2^)(X« + K^Mz) + K^? , (6.27) 

G° = 2p[r^2^2. + rS3^3.]p'(^) + K^p'iz) . (6.28) 

The solution for dX^ can again be written in the form (6.17), but with L^i \ 
/ciand Bi, Ci instead of these expressions with the indice "3". 

Taking into account formulaes (6.25 - 6.28) for Ai, Bi and Ci, the final expression for 
dX^ can be written as 

^p(^)p» V^ip' + + Kf^ + hp" + f2P' + hp + h , , 

dX — , (6.29) 

-^p'{z)JF^p^{z)+F2p{z) + Kf^ +g,p^{z)+g2p{z)+g. 



where Fi, F2, /i, /2, /s, /4, fi-i, 5'2 and g^ are functions, depending on gap, (a, = 1, 2) 
and on the ratios ^, ^, ^, ^, 

Cl ' ai ' (12 ' Cl ' C2 

It is also straightforward to prove that expressions (6.22) for dX'^ and (6. 29) for 
dX^do not represent eUiptic functions. If one assumes that dX^ is an elliptic function, 
then from the standard theory of elliptic functions [9, 11] it will follow that dX^ can be 
represented as 

dX^ = K,{p) + p {z)K2{p) , (6.30) 

where i^i(p) and K2{p) depend on the Weierstrass function only. But this representation 
will contradict with the expression (6.29) for dX^- consequently the initial assumption 
has to be rejected. Similarly, it can be proved that dX'^ is not an elliptic function. The 
details of this simple proof will be left for the interested reader. 

7 COMPLEX COORDINATE DEPENDENCE 
OF THE METRIC TENSOR COMPONENTS 
FROM THE UNIFORMIZATION OF A CU- 
BIC ALGEBRAIC SURFACE 



In this section it will be shown that the solutions (6.17), (6. 22) and (6. 29) of the 
cubic algebraic equation (3.5) enable us to express not only the contravariant metric 
tensor components through the Weierstrass function and its derivatives, but the covariant 
components as well. 

Let us write down for convenience the system of equations (6.17), (6. 22) and (6. 29) 
for dX^, dX'^ and dX^ as (Z = 1, 2, 3) 

dX\X\X\X^)^Fi{g,^{y.),T%{y.),p{z),p{z))^Fi{y.,z) , (7.1) 
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where the appearence of the complex coordinate z is a natural consequence of the uni- 
formization procedure, applied with respect to each one of the cubic equations from the 
"embedded" sequence of equations. 

Yet how the appearence of the additional complex coordinate z on the right-hand side 
of (7.1) can be reconciled with the dependence of the differentials on the left-hand side 
only on the generalized coordinates (X^, X^, X^) (and on the initial coordinates 
because of the mapping X* = X'^ [x^ . x'^ . x'^))l The only reasonable assumption will be 
that the initial coordinates depend also on the complex coordinate, i.e. 

= X\x\z),x'{z),x\z)) = X'(x, z) . (7.2) 



Taking into account the important initial assumptions (Z = 1, 2, 3) 

rIF, 

d^X^ ^O^dFi{X{z),z)^^dz , (7.3) 

dz 

one easily gets the system of three inhomogeneous linear algebraic equations with respect 
to the functions ^ and ^ (/ = 1, 2, 3): 

dFi dX^ dFi dX^ dFi dX^ dFi 

— \ \ \ ^ = {7 A) 

dX^ dz 9X2 dz dX^ dz dz ' ^ ' ' 

The solution of this algebraic system {i,k,l — 1, 2, 3) 

represents a system of three first - order nonlinear differential equations. A solution of 
this system can always be found in the form 

X^^X\z) ■ X2 = X2(z) ; X^^X\z) . (7.6) 

and therefore, the metric tensor components will also depend on the complex coordinate 
z, i.e. Qij = gij{'X{z)). Note that since the functions in the right-hand side of (7.5) 
depend on the Weierstrass function and its derivatives, it might seem natural to write that 
the solution of the above system of nonlinear differential equations gij will also depend 
on the Weierstrass function and its derivatives 

g,j = gij{X\p{z),p{z),X'{p{z),p{z),X'{p{z),p{z)) = gij{z) . (7.7) 

Note however that for the moment we do not have a theorem that the solution of the sys- 
tem (7.5) will also contain the Weierstrass function. But in spite of this, the dependence 
on the complex coordinate z will be retained. 
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8 DISCUSSION 



This paper continues the investigation of cubic algebraic equations in gravity theory, 
initiated in a previous paper [10]. 

Unhke in [10], where the treatment of cubic algebraic equations has been restricted 
only to the choice of the contravariant tensor g'^^ — dX'^dX^ , in the present paper it was 
demonstrated that under a more general choice of g'*-', there is a wide variety of algebraic 
equations of various order, among which an important role play the cubic equations. Their 
derivation is based on two important initial assumptions: 

1. The covariant and contravariant metric components are treated independently, 
which is a natural approach within the framework of affine geometry [15 - 18]. 

2. Under the above assumption, the gravitational Lagrangian (or Ricci tensor) should 
remain the same as in the standard gravitational theory with inverse contravariant metric 
tensor components. 

The proposed approach allows to treat the Einstein's equations as algebraic equations, 
and thus to search for separate classes of solutions for the covariant and contravariant met- 
ric tensor components. It can be supposed also that the existence of such separate classes 
of solutions might have some interesting and unexplored until now physical consequences. 
It has been shown also that the "transition" to the standard Einsteinian theory of grav- 
ity can be performed by investigating the intersection with the corresponding algebraic 
equations. 

The most important result in this paper is given in Section 6 and is related to the possi- 
bility to find the parametrization functions for a multicomponent cubic algebraic surface, 
again by consequent application of the linear-fractional transformation. The parametriza- 
tion functions in this particular case represent complicated irrational expressions of the 
Weierstrass function and its derivative, unlike in the standard two - dimensional case, 
where they are the Weierstrass function itself and its first derivative. The advantage of 
applying the linear- fractional transformations (6.5) and (6.24) is that by adjusting their 
coefficient functions (so that the highest - third degree in the transformation equation will 
vanish), the following sequence of plane cubic algebraic equations is fulfilled (the analogue 
of eq.(65) in [10]): 

P?\H^))m%) + P^'\nm)m%) + (n(3))m(3) + Pg^ = , (8.1) 
(n(2))m^2) + Pf (n(2))m^2) + Pf (n(2))m(2) + Pg^ = , (8.2) 

P['\Hi))ml^ + P2\Hi)Wm + P^Mm^,) + Pg^ = , (8.3) 

where m(3), m(2), m(i) denote the ratios ^, ^, ^ in the corresponding linear - fractional 

~3 - — 2 ~i 

transformations and n(3), n(2), n(i) are the "new" variables dX , dX , dX . The sequence 
of plane cubic algebraic equations (8.1 - 8.3) should be understood as follows: the first one 
(8.1) holds if the second one (8.2) is fulfilled; the second one (8.2) holds if the third one 
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(8.3) is fulfilled. Of course, in the case of n variables (i. e. n component cubic algebraic 
equation) the generalization is straightforward. Further, since each one of the above plane 
cubic curves can be transformed to the algebraic equation {i = 1, 2, 3) 

^(i) = Pi\nii))ml) + P2\n(i))ml) + P3\n(i))m(i) + (n(i)) (8.4) 

and subsequently to its parametrizable form, one obtains the solutions of the initial mul- 
ticomponent cubic algebraic equation. 

Finally, it has been shown that from the expressions (7.1) a system of first - order 
nonlinear differential equations can be obtained, for which always a solution — X^{z), 

= X^{z), X^ = X-^i^z) exists. Thus the dependence on the generalized coordinates 
X^, X^, X"^ in the uniformization functions (7.1 ) dissappears and only the dependence 
on the complex coordinate z remains, as it should be for uniformization functions. 
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